Experimental demonstration of deterministic one-way quantum computing on a NMR 

quantum computer 
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One-way quantum computing is an important and novel approach to quantum computation. By 
exploiting the existing particle-particle interactions, we report the first experimental realization of 
the complete process of deterministic one-way quantum Deutsch-Josza algorithm in NMR, including 
graph state preparation, single-qubit measurements and feed-forward corrections. The findings in 
our experiment may shed light on the future scalable one-way quantum computation. 

PACS numbers: 03.67.Lx, 03.67.Mn, 76.60.-k 



INTRODUCTION 



The one-way quantum computing (QC) [1, 2] is a re- 
cently proposed approach to quantum computation [3]. 
Being entirely different to the traditional quantum cir- 
cuit model [3] , it invokes only single-qubit measurements 
with appropriate feed-forwards to accomplish the com- 
putation, provided a highly entangled state - the cluster 
state [4] or some other special shaped graph states [5] - 
is given in advance. These entangled states serve as the 
universal resource of the one-way QC. The one-way QC is 
not only important as a novel quantum computing model, 
it also helps people to further understand the quantum 
entanglement and measurement since these two funda- 
mental physcial concepts are particularly highlighted in 
this model. 

The one-way QC has so far attracted much attention 
of the physical community. Besides various theoretical 
researches, the existing experiments mainly focused on 
the generation and characterization of a few-qubit graph 
states [6-8], the demonstration of one- and two-qubit 
gates [6, 9-11], and the realization of two-qubit quan- 
tum algorithms [6, 9, 10, 12]. Up to now all the exper- 
iments of one-way QC were performed in linear optics. 
Owing to the lack of interaction between photons, the 
cluster states were generated probabilistically, and the 
success rate decreased exponentially with the number of 
photons. 

In this paper, by exploiting the existing particle- 
particle interactions, we realized the deterministic one- 
way QC in NMR, including the graph state generation, 
single-qubit measurement, and feed-forward. Our exper- 
iment consists of two parts, the deterministic generation 
of a star-like four qubit graph state and the implementa- 
tion of a two-qubit Deutsch-Josza (DJ) algorithm. 
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FIG. 1: (a) Schematic illustration of the Deutsch-Josza algo- 
rithm, (b) The quantum networks for all possible oracles. The 
inter-line symbol in the network denotes the controlled-NOT 
gate |0}(0| (c) <g> I (t) + |l)(l| (c) ® o-x \ where the superscripts 
denote the qubits and I is identity. 



II. THEORY 

We consider the simplest version of the DJ [13] algo- 
rithm, which examines whether an unknown one-bit to 
one-bit function / is constant or balanced. There are 
four possible such functions as described in Fig. 1(b). 
Classically one needs to call the function twice to check 
both the outputs /(0) and /(l), while in DJ algorithm 
only one function call is needed. The process of the al- 
gorithm is illustrated in Fig. 1. The "function call" 
is implemented in the oracle, which applies the follow- 
ing unitary operations: first a <r 2 applied on the tar- 
get qubit (denoted by t), then the two-qubit operation 
\x)c\y)t -> \x) c \y © fj(x)) t (x,y e {0,1}) corresponding 
to the specific function fj (j = 1,2,3,4). The result of 
the algorithm is read on the control qubit (denoted by c) 
by measuring it with a x (the Pauli operator). If the out- 
come is |+) then the function is definitely to be constant, 
otherwise it is balanced. 

To implement the DJ algorithm one should be able to 
construct all the possible four configurations (Fig. 1(b)). 
We find it is sufficient to do these on the star-like 4-qubit 
graph state with appropriate single-qubit measurement 
sequences and corresponding feed-forwards. We start 
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with the introduction of the general logical quantum cir- 
cuit that can be realized on the star-like 4-qubit graph 
(Fig. 2(a)), with arbitrary logical input states and ar- 
bitrary single-qubit measurement bases in the x-y plane. 
For this purpose, as in Ref. [1], we first prepare the 
four physical qubits in the graph into the initial state 
\ipt)i ® |+)2 ® |+)3 <& |'0c)4, where the two arbitrary logi- 
cal input states \ipt) and \i/j c ) of the two logical qubits [6] 
(a target qubit and a control qubit) are initially encoded 
on the physical qubits 1 and 4. Then all physical qubits 
are entangled by the entangling operator 



S = S {12) S {23) 5 (42) , 



(1) 



where = |0)(0|« ® af ] + <g> I (fe) is a 

controllcd-phase gate applied on the physical qubits j 
and k [14]. The logical information is now delocalized. 
It is then manipulated by the measurements carried on 
the physical qubits 1 and 2 with the bases B\{a{) = 
{|ai+)i,|ai_)i} and #2(0:2) = {|"2+)2, |«2— >2> respec- 
tively, where \a±) = ^=(|0)±e lQ |l)) and ai,a 2 are arbi- 
trary angles. After the measurements the logical target 
qubit is transferred to the physical qubit 3, while the log- 
ical control qubit is still encoded on the physical qubit 4. 
The effective logical quantum circuit performed on the 
two logical qubits for the above process is shown in Fig. 
2(b). The si and S2 in the circuit denote the outcome 
of the two measurements, where Sj — 0(1) corresponds 
to the output state \ctj + )(\otj-)) (j = 1,2). The pres- 
ence of the si and s 2 represents the randomness intro- 
duced by the single-qubit projective measurements. In 
one-way QC appropriate feed-forwards [1] compensate 
for this randomness to restore the determinacy of the 
computing. 

The design of the resource entangled state and the 
single-qubit measurement sequences for implementing 
the DJ algorithm is done by comparing the general log- 
ical circuit in Fig. 2(b) with the networks in Fig. 1(b). 
Since the input state in the DJ algorithm is |+) ® |+), 
we set \tpt) = IV'c) = |+)- The resulted entangled state, 
|V>g) = (|+)i|0)2|-)3|+>4+|-)i|l>a|+)3|->4)/\/2, which 
is produced by performing the entangling operator S on 
the initial state |+)i <8 |+) 2 ® |+)3 ® is exactly 

the 4-qubit graph state which corresponds to the star- 
like 4-qubit graph [2] (it is also equivalent to the 4-qubit 
GHZ state under local unitary operations). The mea- 
surement bases and the corresponding feed-forward op- 
erations which are chosen to reproduce the networks of 
all possible oracles are summarized in the Table I. Note 
in the design of the measurements which correspond to 
the constant functions j\ and f 2 , the fact is used that a 
CNOT gate is equivalent to the identity operation when 
it acts on the state |+) <£> |+). The final result of the al- 
gorithm is read by measuring the physical qubit 4. If its 
state (after accounting for the feed-forward operation) is 
|+)(|— )) then the function is constant (blanced). 
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FIG. 2: (a) The star-like 4-qubit graph. The vertices repre- 
sent the four physical qubits, and each edge between them 
indicates a controlled-phase gate applied on the connected 
qubits when to produce the entangled state. The arrows 
marked on vertices 1 and 2 denote the measurements. For 
the detailed description of the initial state, entangling, and 
measurements, see the text, (b) The general logical quantum 
circuit realized on the above graph with arbitrary logical input 
states and arbitrary measurement bases. R n (a) — e~ 
(n = x, z) 
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TABLE I: Measurement bases and feed-forward operations for 
implementing the DJ algorithm. FF^ 3 ' and FF^ 4 ' represent 
the feed-forward operations on physical qubits 3 and 4. 





Measurement bases 


FF (3) 




h 


Si(0),S a (0) 






h 


Bi(0),B 2 (0) 






h 


Bi(tt),B 2 (0) 






h 


Si(7t),S 2 (0) 
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III. EXPERIMENT 

A. Graph state generation 

We used the spins of the four 13 C nuclei in the cro- 
tonic acid dissolved in D2O (see Fig. 3) as the four 
physical qubits, which provide the reduced hamiltonian 

#0 = E^i^/i^ + 2nTi <k Jjkiy ) & ) with the Lar- 
mor angular frequencies u>j and J-coupling strengths 
Jjk- On a Bruker UltraShield 500 spectrometer at 
room temperature, the measured parameters are listed 
in Fig. 3 and their relaxation times are obtained as 
Ti(Cl) = 12.37s, Ti(C2) = 4.89s, T X (C3) = 4.13s, 
Ti(C4) = 4.96s, T 2 (C1) = 376.2ms, T 2 (C2) = 506.7ms, 
T 2 (C3) = 566.5ms, T 2 (C4) = 544.5ms. In experiments, 
we labelled the nuclei spins C2, C4, C3, CI as the physic 
qubits 1, 2, 3, 4 in Fig. 2. 

We first initialized the NMR ensemble to a stan- 
dard pseudopure state with deviation 1 0000) (0000 1 from 
the thermal equilibrium state using the spatial average 
technique [15]. Then a 4-qubit GHZ state \i/j)ghz = 
(|0110) + |1001))/V2 can be created from |0000) by the 
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FIG. 3: Molecular structure of crotonic acid with a table of 
the chemical shifts (on the diagonal) and J-coupling constants 
(below the diagonal). The parameters are given in unit Hz. 
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FIG. 4: Network for implementing the one-way DJ algorithm, 
consisting of the main three stages: (i) graph state generation, 
(ii) single-qubit measurements and (iii) feed-forward correc- 
tion. Here H = e~ l T rTy is the pseudo-Hadamard gate with 
its conjugate H + , and G 2 is the pulsed magnetic filed gradi- 
ent along z axis. In the feed-forward process, the controlled 
gates are performed by the conditionality on the controlling 
spin being in the |0) or |1) state, which are determined by 
the values of (A, B): (A, B) = (0, 1), (1, 1), (1, 0) and (0, 0) 
correspond to the functions f\ to /j. 



network shown in Fig. 4. Here, the controlled-not gate 
CNOT" fe ' could be further decomposed into 
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where i?f } (f ) 



(2) 

-K*/i)<?l k) denotes a § -rotation of the 



qubit k around the y axis, and so forth. The J-coupling 
gate e~ l ~ <J " a * can be realized by the free evolution 
of the spin system under Hq with appropriate rcfocus- 
ing pulses. Those z rotations in CNOT gates can be 
omitted due to the initial state 1 0000) [Hi]. Such a GHZ 
state has the total zero spin quantum number that facil- 
itates us to clean-up the state further by a pulsed mag- 
netic field gradient [17]. Finally, the graph state |-0g) 
is generated from the GHZ state by the local operations 

In experiments, all the singlc-qubit operations are re- 
alized by sequences of strongly modulating NMR pulses 
created by GRAPE algorithm [18]. We maximized the 
gate fidelity of the simulated propagator to the ideal gate, 
taking into account of radio frequency (RF) inhomogenc- 
ity Theoretically the gate fidelities we calculated for ev- 
ery pulse are greater than 0.99, and all the pulse lengths 



FIG. 5: (color online) The ideal (a) and measured (b) density 
matrices of the four-qubit GHZ state \iP)ghz = (|0101) + 
|1010))/a/2. The rows and columns are numbered with the 
computational basis states. 



are 600 [is. 

In order to access the quality of the graph state we 
prepared, we employed the state tomography process 
[19] to obtain the reconstructed density matrix p oxp . 
Since only single-quantum coherences can be directly ob- 
served, the tomography involves 19 repetitions of the ex- 
periment, each with a different readout pulse sequence 
[20] by taking consideration of some unresolved J cou- 
plings. The real part of the reconstructed density matrix 
is shown in Fig. 5 along with the theoretical expecta- 
tion. The imaginary part is small comparable to zero 
in the theoretical expectation. To quantify how close 
the prepared state p exp to the theoretical one pid, we 
used the measure of the attenuated correlation defined 
by c(p exp ) = Tr(p i( ipexp)/ V Tr (Pid)' wn ich takes into 
account both systematic errors and the net loss of mag- 
netization due to random errors and decoherence. The 
value of the correlation for the tomographic readout of 
the GHZ state \i/j)ghz is c(p'ghz) = 0.73, compared to 
the simulated value 0.80 where we adopted the decoher- 
ence model [21] to take the loss of magnetization due to 
spin relaxation. The experimental time for the prepara- 
tion is around 85 ms. To remove the effect of the net 
loss of magnetization due to random errors and spin re- 
laxation, we also calculated the fidelity defined [22] by 

Fghz = Tr(p ld p exp ) I \jTr(p 2 id )Tr(p 2 exp )), which yields 

the value of 0.88. Due to the imperfection of GRAPE 
pulses calculated and the effect of strongly-coupling, the 
total theoretical fidelity is around 0.92. Moreover, the 
lower experimental fidelity of 0.88 is mainly caused by 
other uncertainties (e.g., the imperfection of the static 
magnetic field) in the experiments. The correlation of 
the graph state \iPg) is similar to the state \4>Ghz), be- 
cause it was created only by local operations. 

Moreover, we confirmed the generation of the four- 
partite pseudo-entanglement by obtaining a negative ex- 
pectation value of the entanglement witness [23] W = 
1/2 — |V'Ghz)('0Ghz| (I denotes identity operator) from 
the tomographic readout: Tr(VFp C xp) = —0.175. 
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B. Implementation of DJ algorithm 

The required single-qubit projective measurements in 
the one-way QC are absent in current NMR quantum 
computing technologies. But instead wc can use the 
pulsed magnetic field gradients to mimic them [24] , whose 
effects are equivalent to that by applying these projective 
measurements on every member of the ensemble. For ex- 
ample, one could use the operations 

pW . Gl _ #(2)^) _ jjWfr) _ d - R^(ir) - i?W(7r) 

-Gx - flL a J(7r) - P^(tt) - Gi - P^(tt) - R^li-K) 

(3) 

to mimic the measurement on the physical qubit 
1, where Gi is the pulsed magnetic field gradient along 
z with a period of ri/4. These operations dephase all 
the coherences which are associated with the transitions 
of qubit 1 in the density matrix p = | < &)( ( 3 ) | (suppose 
1$) = |O)i|0o) + |l)i|<M, where |</> (i)) is the state of 
the other three qubits), resulting in an ensemble-average 
density matrix |O)i(O||(fo)(0o| + |l)i(l||0i)(</>i|, which is 
exactly the same as the outcome by performing projec- 
tive measurement on every member of the ensemble. 
To mimic the measurement a4 , one just first rotate the 
qubit 1 to the z axis before Pi . In our experiment 

we use the sequences P^(f ) - P z (1) and P^d ) - 
to mimic the measurements with the bases Pi(0) and 
B\ (jr) respectively. Therefore the different outcomes are 
labeled in the different subspaccs of qubit 1 denoted by 
|0)i(0|(i.e., si = 0) and |l)i(l| (i.e., si = 1). The mea- 
surement on qubit 2 with basis Pi(0) is mimicked by a 
similar sequence R_y(^) — Pz 2 \ where 

Pj 2 ) : if )(tt) ~ R { y 4) M - G 2 - P«(7r) - PW(tt) - G 2 

-P^« - R { _%) - G 2 - P«(tt) - flW(Tr) - G 2 , 

(4) 

where G 2 is the pulsed magnetic field gradient with a 
period of t 2 /4. 

After the measurements, the computational outcomes 
have been stored in the different subspaces depending 
on the measurement results si and s 2 . However, to ob- 
tain a deterministic and correct computation, the feed- 
forward operations should be carried out conditionally on 
the measurement results si and s 2 . We also experimen- 
tally realized the feed-forward operations by the condi- 
tional unitary operations illustrated in Fig. 4. After this, 
the judgement of whether the function / is constant or 
balanced is determined by measuring the state of qubit 4: 
if the state of qubits is in |+), the function / is constant; 
conversely, if the state of qubits is in |— ), the function / 
is balanced. As a result of |±) (±| = ^I±a x , the |+)/|— ) 
state will give a positive/negative NMR signal if we first 
reference the spectrum of the thermal equilibrium after 
a [j] y pulse as a positive one. The results are shown in 




"Z. 80 40 -40 -80 ~80 40 : 4(j ^S0 



NMR Frequency (Hz) 

FIG. 6: (color online) The experimental (line) and simulated 
(dotted line) spectra of qubit 4 for one-way DJ algorithm: 
(a)-(b) for the constant functions /i and /2, (c)-(d) for the 
balanced functions fz and fi. 



Fig. 6 for four cases f\ — f±. In (a) and (b), the signal 
amplitude is positive, which indicates the state of qubit 
4 is in |+) and the function / is constant. In (c) and (d), 
the signal is inverted which indicates the state of qubit 
4 is in |— ) and the function / is balanced. Here the 
carbon signal-to-noise ratio (Signal/Noise) of the spectra 
is about 44.4-48.3. The imperfection of the experimen- 
tal spectra compared to the simulated spectra is mainly 
caused by the unideal input graph state, the imperfection 
of the static magnetic field and the GRAPE pulses. Al- 
ternatively, the feed-forward process can also be replaced 
by the single-qubit measurement under a suitable basis 
[1], which makes that the one-way QC can be completed 
solely by single-qubit measurements. 



IV. CONCLUSION 

We realized the complete process of deterministic one- 
way QC using a star-like 4-qubit graph state in NMR: 
the graph state was generated deterministically by ex- 
ploiting the existing particle-particle interactions in our 
experiments; the single qubit measurements were mim- 
icked by pulsed magnetic field gradients; and the feed- 
forward corrections were realized by conditional-unitary 
gates in our specific NMR system. Here NMR technique 
is used as a test-bed for the demonstration of determin- 
istic one-way QC, however, it does not imply that any 
of the results obtained are specific to NMR. The expe- 
rience obtained from our experiments could be helpful 
to other physical systems for future scalable determin- 
istic one-way QC and deserve further investigation. If 
the direct particle-particle interactions of the quantum 
system (such as trapped ions and optical lattices) are 
provided, one can prepare the cluster state determinis- 
tically, do single-qubit measurements, and perform the 
feed-forward correction by unitary gates, therefore real- 
izing scalable one-way QC deterministically. 
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